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Abstract 

The  method  of  state-space  averaging  has  been  successfully 
applied  to  pulse-width  modulated  power  converters,  but  has 
its  limitations  with  switched  circuits  that  do  not  satisfy  a 
“small  ripple"  condition.  This  work  considers  a  more  general 
averaging  procedure  that  encompasses  state-space  averaging 
and  is  potentially  applicable  to  a  much  broader  class  of  cir¬ 
cuits  and  systems.  In  particular,  the  technique  is  shown  to 
be  effective  on  a  number  of  examples  including  resonant  type 
converters. 

1  Introduction 

State-space  averaging  [1,  2,  3]  has  been  demonstrated  to  be 
an  effective  method  for  analysis  and  control  design  in  pulse- 
width  modulated  (PWM)  switching  power  converters.  How¬ 
ever.  as  has  been  noted  in  the  literature,  the  class  of  convert¬ 
ers  that  this  method  can  be  applied  to  is  limited.  Conditions 
for  the  justification  of  state-space  averaging  have  been  char¬ 
acterized  by  a  “small  ripple”  condition,  by  a  “linear  ripple” 
approximation  [1],  and  by  the  degree  to  which  certain  vec¬ 
tor  fields  commute  [2].  With  the  small  ripple  approximation, 
the  assumption  is  that  a  Fourier  series  expansion  for  a  finite 
length  segment  of  a  circuit  waveform  should  be  dominated 
by  its  DC  term.  The  linear  ripple  approximation  requires 
that  the  circuit  waveforms  appear  to  be  linear  functions  of 
time  when  examined  over  a  time  interval  inbetween  switch 
instances.  This  condition  has  been  stated  more  precisely  in 
(2)  in  terms  of  the  Lie  bracket  of  a  pair  of  matrices.  A  re¬ 
cent  paper  [6]  studied  the  application  of  classical  averaging 
techniques  in  PWM  circuits,  and  has  developed  an  asymp¬ 
totic  framework  where  solutions  to  switched  systems  may  be 
approximated  to  arbitrary  accuracy  by  a  power  series  in  a 
small  parameter.  The  small  parameter  is  evidently  related  to 
the  ratio  of  the  switching  period  and  the  system  time  con¬ 
stants.  We  note  here  that  this  ratio  is  typically  small  for  fast 
switching  PWM  circuits,  but  this  is  typically  not  the  case  for 
resonant  type  converters. 

Because  of  the  conditions  listed  above  that  limit  the  ap¬ 
plicability  of  state-space  averaging,  the  method  cannot  be 
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applied  to  a  wide  range  of  power  circuits  that  includes  the 
resonant  type  converters.  The  basic  limitation  in  resonant 
converters  is  that  these  circuits  have  state  variables  that  ex¬ 
hibit  predominantly  oscillatory  behavior.  This  paper  investi¬ 
gates  a  more  general  averaging  scheme  that  can,  in  principle, 
accomodate  arbitrary  types  of  waveforms.  The  method  is 
based  on  a  time-dependent  Fourier  series  representation  for  a 
“sliding  window"  of  a  given  waveform.  For  example,  for  an 
arbitrary  time-domain  waveform  x(»),  the  method  considers 
the  Fourier  coefficients  of  x(s)  for  a  €  (t  —  T,  1]  at  the  time 
instant  t.  Simplifying  approximations  can  be  made  by  omit¬ 
ting  insignificant  terms  in  this  series.  For  instance,  to  recover 
the  traditional  state-space  averaged  model,  one  would  retain 
only  the  DC  coefficient  in  this  averaging  scheme. 

Previously  developed  methods  for  analysis  and  control  de¬ 
sign  in  resonant  converters  have  relied  mainly  on  two  ap¬ 
proaches,  sampled-data  modeling  [4,  5]  and  phase-plane  tech¬ 
niques  [7,  8).  The  sampled-data  approach  taken  in  [4,  5]  re¬ 
sults  in  a  small-signal  model  for  the  underlying  resonant  con¬ 
verter  with  the  perturbation  in  switching  frequency  as  the 
input.  One  difficulty  with  this  approach  is  the  requirement 
of  obtaining  a  nominal  periodic  solution  as  a  first  step  in  the 
analysis.  The  utility  of  the  resulting  model  is  then  limited  by 
its  small  signal  nature.  The  phase  plane  method  of  [7]  is  a 
basic  approach  to  obtaining  a  steady  state  solution  for  a  res¬ 
onant  converter,  and  the  control  scheme  of  [8J  based  on  this 
method  is  evidently  effective.  A  limitation  of  this  approach 
is  its  restriction  to  second  order  systems;  it  is  not  obvious 
how  one  can  incorporate  additional  state  variables  that  are 
associated  with  the  load  or  the  source  dynamics.  In  contrast, 
the  method  presented  in  this  paper  is  an  averaging  scheme 
that  can  be  applied,  in  principle,  to  any  periodically  driven 
system. 

The  remainder  of  the  paper  is  organized  as  follows.  Sec¬ 
tion  2  reviews  some  of  the  main  techniques  used  in  analyzing 
oscillations  in  nonlinear  systems.  Our  approach  is  outlined 
in  Section  3  and  will  be  seen  to  be  very  closely  related  to 
the  classical  methods  discussed  in  Section  2.  Examples  of  the 
application  of  our  method  in  resonant  and  in  PWM  convert¬ 
ers  are  presented  in  Section  4.  Appendix  A  contains  some 
comments  on  the  justification  of  our  scheme. 

2  Averaging  and  Oscillations  in 
Nonlinear  Systems 

One  of  the  pioneers  in  the  study  of  oscillalions  in  nonlin¬ 
ear  systems  was  van  der  Pol  [11]  who  studied  a  periodically 
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driven,  oonlinearly  damped  oscillator  modeled  by 

t>  +  <(>(v)v  +  =  u>lBstn(u)\t)  (1) 

The  paper  [11]  postulated  a  solution  of  the  form 

v(t )  =  b,(t)sin(u/tt)  +  />j(t)cos(u>it)  (2) 

where  b^(t)  and  b^{t)  were  supposed  to  be  slowly  varying  func¬ 
tions  of  time.  That  is,  for  /b  =  1,2, 

M<)  <  (3) 

i*( 0  <  wih»(0  (<) 

Van  der  Pol’s  method  proceeds  by  firstly  inserting  the  pos¬ 
tulated  solution  (2)  into  (1)  and  extracting  the  fundamental 
harmonic  response  from  the  nonlinear  element  in  the  model 
( 1).  (The  extraction  of  the  fundamental  harmonic  response  is 
known  as  harmonic  linearization  and  is  essentially  equivalent 
to  the  describing  function  method  described  below.)  Sec¬ 

ondly.  harmonic  balance  is  imposed  to  obtain  two  equations, 
one  in  sjn(u-'il)  and  one  in  eos(a;|f).  The  result  is  a  state-space 
model  that  governs  the  behavior  of  the  coefficients  &i(<)  and 
The  method  *o  be  introduced  in  Section  3  draws  on 
these  ideas,  but  makes  slightly  different  assumptions. 

A  similar  approach  to  the  study  of  the  second  order  sys¬ 
tem  (1)  was  taken  in  the  work  of  Krylov,  Bogoliubov,  and 
Mitropolsky  [12].  This  work  relied  on  the  assumption  of 
slowly  varying  phase  and  amplitude,  i.e.  postulated  a  so¬ 
lution  of  the  form 

t>(f)  =  a(t)sin{uJot  +  0(f)}  (5) 

t>(<)  =  a(t)wocos{u.-0l  +  0(0}  (6) 

where  a  1)  and  8(t)  were  supposed  to  be  slowly  varing  func¬ 
tions  of  time.  An  implicit  constraint  arising  from  (5,6)  IS 
that 

a(t).«m{u.0<  +  #(<)}  +  n(<)0(<)cos{wo<  +  0(0}  =  0-  (") 

The  method  of  [12]  proceeds  similarly  to  van  der  Pol’s  tech¬ 
nique  by  substituting  the  postulated  form  of  solution  (5,6) 
into  the  nonlinear  differential  equation  (1)  to  derive  a  state- 
space  equation  for  the  variables  a(t)  and  0(0-  The  assump¬ 
tion  of  slowly  varying  phase  and  amplitude  is  used  to  deduce 
an  or traged  model  for  the  evolution  of  a(t)  and  0(1).  The  av¬ 
eraging  operation  performed  in  the  method  of  [12J  is  actually 
equivalent  to  harmonic  linearization  in  that  the  fundamental 
component  of  the  response  of  the  nonlinear  element  of  (1)  is 
extracted.  See  [12]  for  more  details  on  this  method  or  the 
summaries  in  [10,  13]- 

Anotber  technique  that  is  closely  related  to  the  methods 
outlined  above  is  the  describing  function  method  discussed 
in  detail  in  Gelb  and  Vander  Velde  [13].  This  method  relies 
on  harmonic  balance  for  analyzing  steady  state  and  transient 
oscillations  in  nonlinear  systems.  To  understand  the  method, 
consider  the  autonomous  system  of  Figure  1  that  has  a  static 
nonlinearity  N(»)  in  the  feedback  path.  The  input  to  the 
nonlinearity  is  labeled  r(i)  and  the  output  is 

y(t)  =*{*(!)}.  (8) 

The  sinusoidal  describing  function  method  is  a  technique  for 
analyzing  limit  cycles  in  a  system  of  the  type  of  Figure  1  or 


Figure  1:  Nonlinear  Feedback  System 

in  such  a  system  that  has  an  external  periodic  input.  The 
basic  premise  of  the  describing  function  method  is  that  r(() 
is  approximately  a  sinusoidal  waveform,  i.e. 

x(t)  *3  ,Ystn(mi).  (9) 

If  the  nonlinearity  is  not  too  severe,  that  is  it  does  not  gen¬ 
erate  excessively  large  harmonics,  and  if  the  linear  transfer 
function  G( s)  is  a  good  low  pass  filter,  the  assumption  that 
x (t)  is  sinusoidal  is  quite  plausible.  The  describing  function 
method  characterizes  y(t )  by  its  fundamental  frequency  com¬ 
ponent 

»{t)*/f(X)x(t)  (10) 

where  A'{X)  is  the  ratio  of  the  fundamental  component  of 
y[t)  to  that  of  x(f).  -V(.Y)  is  called  the  sinusoidal  describing 
function  for  Ar(»)  and  is  dependent  on  the  amplitude  X  of 
the  oscillation.  For  the  loop  of  Figure  1,  the  condition  for  a 
steady  state  oscillation  is 

1  +  GO'w)V(A')  =  0,  (11) 

or  equivalently 

G(ju;)  =  -1/A'(.Y).  (12) 

The  power  of  the  method  is  evident  from  (12)  since  this 
equation  can  be  solved  graphically  by  plotting  G(ju>)  and 
—  lfAf(X)  on  the  complex  plane  and  determining  intersec¬ 
tion  points  of  these  graphs.  Each  intersection  point  (.Y*,u»*) 
generates  a  candidate  for  an  approximate  limit  cycle  solution 
of  the  form 

*(f) »  X-sin(v't).  (13) 

The  paper  of  Bergen  and  Franks  [14]  gives  conditions  under 
which  an  exact  periodic  solution  is  nearby  and  gives  error 
bounds  on  the  approximate  solution  in  terms  of  the  opera¬ 
tors  G(jw)  and  /V(e).  The  describing  function  method  can 
also  be  used  when  periodic  external  signals  enter  the  system. 
In  this  case,  one  seeks  a  periodic  solution  with  fundamental 
frequency  equal  to  the  drive  frequency.  The  sinusoidal  de¬ 
scribing  function  method  generates  the  amplitude  and  phase 
of  an  approximate  solution. 

The  describing  function  technique  has  been  generalized  for 
the  analysis  of  transient  oscillations  in  [13].  The  generaliza¬ 
tion  assumes  a  solution  of  the  form 

x(<)  =  .Y(l)e*M,  (14) 
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which  cm  be  identified  with  the  method  of  slowly  varying 
phase  and  amplitude  (12)  if  one  takes  4>(t)  =  +  0(t).  The 

transient  analysis  proceeds  by  noting  that 

i(t)  =  -YeJ*<,)  + 

=  (Y  +  (15) 

and  identifying  these  quantities  by 

i(t)  =  ( <r(t )  +  jw(t))i(t)  =  s*(t).  (16) 

Although  not  clearly  stated  in  (13),  the  method  evidently 
then  replaces  G(jui)  by  G(ff  +  ju)  and  imposes  harmonic 
balance  by  replacing  nonlinearities  with  amplitude-dependent 
describing  functions.  The  result  is  a  pair  of  coupled  differ¬ 
ential  equations  (corresponding  to  real  and  imaginary  com¬ 
ponents)  in  <r(t),uj(t),  and  the  derivatives  of  these  variables. 
Reference  (13)  continues  by  giving  approximate  schemes  for 
solving  the  obtained  differential  equations.  Our  method  con¬ 
sidered  in  the  sequel  is  closely  related  to  the  transient  and 
steady  state  analysis  techniques  that  rely  on  describing  func¬ 
tions. 

3  Generalized  Averaging 

Our  averaging  method  is  based  on  the  fact  that  the  waveform 
r(*)  can  be  approximated  on  the  interval  (t  —  T,  <]  to  arbitrary 
accuracy  with  a  Fourier  series  representation  of  the  form 

x(f. »)  =  x(t  -  T  +  s)  =  £»4(t)c-’*“*,‘-T+‘>  (17) 

* 

where  the  sum  is  over  all  integers  k,  w,  =  2x/T,  3  €  (0,T), 
and  the  (-r)k(O  are  complex  Fourier  coefficients.  These 
Fourier  coefficients  are  functions  of  time  and  are  determined 

by  r 

(*)*(« )  =  f  l  -  T  +  3)e-jfc"-(,-r+*>  <L 3.  (18) 

The  k-th  coefficient  will  also  be  referred  to  as  the  index-k  co¬ 
efficient.  Note  that  the  strategy  here  is  to  consider  a  window 
of  length  T  of  a  given  waveform,  and  to  view  the  waveform 
as  periodic  (or  having  finite  duration  T)  so  that  it  can  be 
analyzed  with  a  Fourier  series.  The  remaining  steps  in  the 
analysis  simply  compute  the  time-evolution  of  these  Fourier 
series  coefficients  as  the  window  of  length  T  slides  over  the 
actual  waveform.  Our  approach  in  doing  so  is  to  determine  an 
appropriate  state-space  model  in  which  the  coefficients  (18) 
are  the  states.  In  this  way,  the  technique  is  very  similar  to 
the  classical  methods  for  analyzing  nonlinear  oscillations. 

As  an  aside,  note  that  one  possible  approach  for  deriving 
the  theory  of  state-space  averaging  is  to  consider  the  one-cycle 
average 

At)  =  jj'_T*(3)d3  (19) 

for  the  state  x(<)  of  a  switching  converter  operating  at  fre¬ 
quency  1  IT.  The  connection  with  our  scheme  is  that  7(f)  — 
(•r)o(f)  corresponds  to  the  DC  coefficient  in  the  Fourier  series 
representation  (17).  More  details  on  this  are  given  Section  4. 

Certain  properties  of  the  Fourier  coefficients  (18)  are  key 
for  the  analysis  and  are  detailed  l>elow. 


Differentiation  with  respect  to  time  The  time  deriva¬ 
tive  of  the  k-th  coefficient  is  computed  to  be 

(/)  -  jku,(x)k{t).  (20) 

This  formula  will  be  very  important  in  computing  the  form 
of  an  averaged  model  involving  the  Fourier  coefficients.  The 
case  where  u>,  is  time  varying  will  also  need  to  be  considered 
for  the  Malysis  of  systems  where  the  drive  frequency  is  not 
constant,  e.g.  in  resonant  type  converters.  In  this  case,  the 
formula  (20)  is  only  m  approximation,  but,  for  slowly  varying 
wa(t),  it  is  a  good  approximation.  Appendix  A  gives  details 
on  the  approximations  involved. 

Transforms  of  functions  of  variables  Another  impor¬ 
tant  ingredient  in  the  Fourier  series  coefficient  representation 
of  a  signal  involves  the  computation  of  the  following: 

{/(xi,x3,...,r,))4  (21) 

where  /(•,...,•)  is  a  general  scalar  function  of  its  arguments. 
One  approach  for  computing  the  quantity  (21)  is  by  the  de¬ 
scribing  function  method  (13).  As  outlined  earlier,  this  would 
require  the  assumption  that  the  arguments  X|,  are 

sinusoidal  in  nature  or  perhaps  have  some  other  fixed  wave¬ 
shape.  For  instance,  one  special  case  is  where  it  is  desired  to 
compute  (/(x)),  and  (/(xj)_i,  and  it  is  known  that  i  is  dom¬ 
inated  by  its  fundamental  Fourier  components.  This  would  fit 
the  framework  of  the  sinusoidal  describing  function  method. 
An  example  of  this  type  is  given  in  Section  4  in  the  context 
of  a  resonant  converter. 

A  procedure  for  exactly  computing  (21)  is  available  in  the 
case  where  /(•,  ...,*)  is  polynomial.  The  procedure  is  based 
on  the  following  convolutional  relationship: 

(xy)t  =  £(r)*-i(v><  (22) 

i 

where  the  sum  is  taken  over  »’]  integers  «.  In  many  cases,  we 
shall  rely  on  many  of  the  terms  in  the  series  in  (22)  being  neg¬ 
ligibly  small.  The  quantity  in  (21)  can  be  computed  in  the 
case  where  /(#,...,•)  is  a  polynomial  by  considering  each 
homogeneous  term  separately.  The  constant  and  linear  terms 
'-r  rivial  to  transform.  The  transforms  of  the  quadratic 
uni'  can  be  computed  using  (22).  Homogeneous  terms  of 
h.  '  r  jrder  cm  be  dealt  with  by  factoring  each  such  term 
in*.  e  product  of  two  lower  order  terms.  Then  the  proce¬ 
dure  cm  be  applied  to  each  of  the  factors.  This  process  is 
guaranteed  to  terminate  since  factors  with  only  linear  terms 
will  eventually  arise. 

Application  to  state-space  models  of  power  electronic 
circuits  Here  we  consider  the  application  of  the  method  to 
a  state-space  model  that  has  some  periodic  time-dependence. 
This  is  the  type  of  model  that  typically  arises  in  a  switched 
power  electronic  circuit.  For  instMce,  consider  the  model 

|x(t)  =  /{x(t),u(0)  *  (23) 

where  u(t)  is  some  periodic  function  of  time  with  period  T. 
The  variable  u(t)  may  be  the  ramp  function  used  to  imple¬ 
ment  a  PWM  scheme,  or  may  be  a  square  wave  of  source  volt¬ 
age  applied  to  a  resonant  tank  circuit  in  a  resonant  converter. 
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To  apply  the  generalized  averaging  scheme  to  a  converter  with 
model  (23),  we  simply  compute  the  relevant  Fourier  coeffi¬ 
cients  of  both  sides  of  (23),  i.e. 


L  =  197  pH 


C  =  lOOnF 


{ltX)t  =  m  X’U)>‘  (24) 

for  the  k- th  coefficient.  A  first  step  in  simplifying  the  model 
is  to  incorporate  the  rule  for  computing  the  derivative  of  the 
fc-th  coefficient.  We  obtain 

|(x)*  = ->*«.(*)*  +  (/(*,  u))».  (25) 

The  second  term  on  the  right-hand  side  of  (25)  can  be  sim¬ 
plified  into  explicit  functions  of  the  coefficients  (x),  and  (u}< 
using  describing  functions.  The  essence  in  modeling  is  to  re¬ 
tain  only  the  relatively  large  Fourier  coefficients  to  capture 
the  interesting  behavior  of  the  system.  As  previously  dis¬ 
cussed.  we  would  retain  only  the  index-zero  (DC)  coefficients 
for  a  fast  switching  PWM  circuit  to  capture  the  low  frequency 
behavior.  The  result  would  be  precisely  the  state-space  aver¬ 
aged  model.  For  a  resonant  converter  which  has  some  of  its 
states  exhibiting  predominantly  fundamental  frequency,  sinu¬ 
soidal  waveforms  and  some  of  its  states  exhibiting  predom¬ 
inantly  DC  (or  slowly  varying)  waveforms,  we  would  retain 
the  index-one  (and  minus  one)  coefficients  for  those  states 
exhibiting  sinusoidal-type  behavior  and  the  index-zero  coeffi¬ 
cients  for  those  states  exhibiting  slowly  varying  behavior. 

The  method  is  demonstrated  on  some  examples  in  the  fol¬ 
lowing  section. 

4  Examples 
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Series  Resonant  Converter  with  Voltage  Source  Load 
The  first  example  to  be  considered  is  the  series  resonant  con¬ 
verter  of  Figure  2  that  has  a  voltage  source  load.  The  circuit 
parameters  are  as  given  in  the  figure  and  the  diodes  are  as¬ 
sumed  to  be  ideal.  Note  that  the  resonant  tank  frequency  is 
approximately  36KHz. 

A  state-space  model  for  this  circuit  takes  the  form 


-{-v  -  V„3gn(i)  +  V'<3sn(sin(w,l))} 


C‘ 


(26) 


where  syn(s)  indicates  the  sign  function.  Typical  waveforms 
for  this  circuit  are  shown  in  the  figure.  These  waveforms  were 
generated  by  stepping  the  drive  frequency  between  38KHz 
and  tOKHz.  One  of  the  important  features  of  the  waveforms 
is  that  following  each  step  change  in  the  driving  frequency, 
the  waveforms  appear  to  be  amplitude  modulated  sinusoids. 
The  waveforms  evidently  settle  down  to  an  approximately  si¬ 
nusoidal  steady  state.  For  this  reason,  these  waveforms  may 
be  very  well  approximated  with  the  fundamental  frequency 
term;  in  the  Fourier  series  (17).  We  are  led  to  examine  a 
model  containing  only  the  coefficients  (i)i,  (»)-i,  (v)i,  and 
(r)_i .  This  can  be  obtained  by  considering  the  application  of 
the  operator  (*)i  (or  equivalently  (*)-i)  to  the  model  (26), 
i.e. 


Figure  2:  DC-DC  Series  Resonant  Converter  with  Voltage 
Source  Load 


^(«),  =  -jw.(«),  +  j-{-(o),  -  V,{sgn{i))i  + 
V,(sgn(sin(u>,t))),} 

jj(v)i  =  -MW,  +  ^(i)i  (27) 

Note  that  the  elements  of  the  two  component  state  vector  of 
(27)  are  complex  Fourier  coefficients,  and  so  this  model  actu¬ 
ally  corresponds  to  a  fourth-order  state-space  model  with  real 
variables.  The  corresponding  real  fourth-order  model  could 
be  obtained  equivalently  from  the  index-(-l)  coefficients.  If 
the  window  length  is  taken  to  be  T  =  2 */u>„  the  term 

(sgn(sin(ui,i)))x  =  (28) 

is  simply  the  constant  amplitude  of  the  first  coefficient  of  the 
Fourier  series  for  a  square  wave.  The  term  (-."onfi)),  can  be 
evaluated  using  the  describing  fi-iict’ou  appioach  by  assuming 
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t(( )  is  approximated  as  a  sinusoid  over  each  interval  of  length 
T.  In  this  case,  we  have 

<s,n(,)),  =  V<'>*.  (29) 

Consequently,  the  model  (27)  is  approximated  with  the  time- 
invariant  model 

£(«).  =  -MW, +  i{-< v)x-Ve\^  -V.jl) 

|(0,  =  -jw.(u),  +  !(.),.  (30) 

The  variables  u>,  and  V,  may  be  considered  as  external  in¬ 
puts  to  this  model.  One  question  that  arises  is  how  well  the 
magnitudes  of  the  complex  coefficients 

I  (Oil  =  v/(^(‘).),  +  (/m(.),)i 

||{t»),||  ='v/(^(,-)1)1  +  (/mW.)i.  (31) 

coincide  with  the  amplitude  of  the  waveforms  in  Figure  2. 
The  simulation  of  Figure  3  compares  precisely  these  quanti- 


to.  i  and  ||(i)l|| 


0.  7. SI-4  0.0015  0.00J2S  ^0.001 

Figure  3:  Amplitude  of  First  Coefficients  Compared  to  Actual 
Waveforms 

ties.  This  figure  plots  on  the  same  axes  the  waveforms  v(f) 
and  ||(r)i(<)||  for  the  case  where  u>,  undergoes  a  step  cha..ge, 
and  the  figure  similarly  plots  i(f)  and  ||(i)i(()||-  It  is  evident 
from  tbe  figure  that  the  correlation  between  the  two  types  of 
waveforms  is  excellent. 

The  utility  of  the  method  can  now  be  realized.  For  instance, 
it  is  straightforward  to  obtain  a  steady  state  solution  for  this 
model  by  setting  the  derivatives  of  (u)i  and  («)i  to  zero.  This 
solution  can  be  computed  in  a  closed  form.  After  obtaining 
a  steady  state  solution,  the  model  may  be  linearized  about 
the  steady  state  to  obtain  small  signal  transfer  functions  from 
inputs  such  as  switching  frequency  w,  or  source  voltage  V,  to 
variables  such  as  (v)i  or  (i)j.  For  purposes  of  illustration,  we 
carry  out  this  procedure  here. 

The  steady  state  solution  is  determined  to  be 


By  considering  perturbations  around  this  steady  state  solu¬ 
tion  it  is  possible  to  compute  the  transfer  function  from  in¬ 
put  switching  frequency  w,  to  the  capacitor  voltage  amplitude 
il(u)i||.  This  is  given  by 

w  A'a  +  2A 

~Mu'V  +  Ar»  +  2EsJ  +  AEs  +  A1  (34) 

where  M  =  |<»)rl,  K  =  2u>lV./wMu>„  =  1  /LC,  £  = 
Uo  +  and  A  =  uj  —  w*.  A  root  locus  of  the  poles  of  this 
transfer  function  is  shown  in  Figure  4.  The  plot  shows  the 
location  of  the  poles  as  tbe  load  voltage  V„  increases  from 
zero  to  V,.  The  poles  are  seen  to  be  damped  for  0  <  Vt  < 


Figure  4:  Root  Locus  of  Small  Signal  Poles  as  Load  Voltage 
Varies  from  Zero  to  V, 

V,.  For  small  values  of  V„  there  is  a  pair  of  poles  located 
near  ± j (u>,  —  u>o).  This  pair  evidently  governs  the  relatively 
slow  amplitude  modulation  seen  in  Figure  3.  One  argument 
for  the  presence  of  this  slow  amplitude  modulation  is  that 
it  is  the  result  of  beating  between  oscillations  at  the  natural 
frequency  un>  and  the  driving  frequency  u>,.  Note  that  the 
locus  of  Figure  4  closely  resembles  the  locus  of  reference  [4] 
that  was  determined  using  an  “exact"  sampled-data  model. 
It  is  also  of  interest  that  the  transfer  function  (34)  has  a 
right-half  plane  zero  for  A  <  0,  that  is  for  frequencies  below 
resonance. 

One  interesting  feature  of  the  simulation  of  Figure  3  is  that 
the  high  frequency  modes  of  the  model  (34)  are  not  evident  in 
the  transient  response.  This  can  be  attributed  to  the  fact  that 
these  modes  are  weakly  coupled  to  the  input  corresponding 
to  perturbations  in  driving  frequency  u>t.  • 

DC-DC  Series  Resonant  Converter  with  Capacitor 
Load  A  more  complex  example  is  the  series  resonant  con¬ 
verter  with  a  capacitor  load  that  is  shown  in  Figure  5.  For 
this  example,  a  state-space  model  can  be  written  as 

d  \ 

7l'  =  £ ;{-t’-*VJn(0  + V,spn(3rn(w,0)} 

d  I  . 


Figure  5:  DC-DC  Series  Resonant  Converter  with  Capacitor 
Load 

Here.  «£-;(• )  and  a b»(»)  are,  respectively,  the  sign  and  abso¬ 
lute  value  functions. 

In  our  analysis,  we  retain  the  fundamental  frequency 
Fourier  coefficients  for  i  and  v  and  the  DC  coefficient  for 
r0.  This  choice  can  he  motivated  by  studying  the  converter 
waveforms  and  noting  that  the  variables  i  and  v  exhibit  pre¬ 
dominate  ly  sinusoidal  behavior  while  the  output  voltage  v0 
exhibits  rather  slowly  varying  dynamical  behavior.  If  this  was 
not  the  case,  more  terms  could  be  retained  in  the  averaged 
model.  The  key  steps  in  applying  the  averaging  operation 
to  the  model  (3-5)  are  in  computing  {tgn(i)v„)i  and  (a6s(i))0. 
These  can  be  determined  using  describing  function  methods  if 
i  and  u  are  considered  to  be  sinusoidal  and  v0  is  approximated 
with  a  constant  term.  With  these  simplifying  assumptions, 
we  find 


( *9" (* K)  1  85 

(36) 

(abs(i))o  »  ~||(*>i H- 

(37) 

In  order  to  verify  the  obtained  average  model,  a  simulation 
that  compares  waveforms  generated  by  this  model  with  wave¬ 
forms  of  the  underlying  system  is  shown  in  Figure  6.  As 
is  evident  from  the  figure,  the  averaged  model  predicts  the 
transient  behavior  quite  accurately.  Analysis  steps  similar  to 
those  carried  out  in  the  previous  example  could  be  performed, 
but  we  omit  the  details.  • 


PWM  Up-Down  Converter  In  this  final  example,  we  il¬ 
lustrate  how  to  refine  the  method  of  state-space  averaging  by 
including  higher  order  terms  in  the  Fourier  series  expansion 
(17).  The  converter  of  Figure  7  operating  in  continuous  con¬ 
duction  mode  can  be  modeled  with  a  second  order  state- space 
model  of  the  form 

~t  =  Ai  +  uBz  +  bu  +  f  (38) 


0.  0.001)  0.0022)  0.00) 

Figure  6:  Comparison  of  Transient  Waveforms  of  Avaraged 
Model  with  Waveforms  of  Underlying  Resonant  Converter 
with  Capacitor  Load 


/,  =  10A Hz 
£>  =  0.5 


Figure  7:  Up-Down  Converter 


where  the  variable  u  takes  the  values  0  and  1  which  corre¬ 
spond  to  the  instaneous  switch  position.  For  nominal  PWM 
operation  at  the  fixed  frequency  /,  =  u>,/2 x,  the  state-space 
averaged  model  can  be  obtained  by  applying  the  one-cycle  av¬ 
eraging  operation  (*)o  to  the  model  (38).  A  refinement  of  the 
slate-space  averaged  model  can  be  obtained  by  considering 
an  additional  term  in  the  series  (17)  that  corresponds  to  the 
fundamental  component  of  the  ripple.  For  this,  we  consider 
the  application  of  both  (•)<>  and  (*)i  to  the  model  (38).  The 
result  is 

jt(x)0  =  A{x)  o  +  P((x)o(u)0  +  (*)i(«)-i  +  (-r)-i(u)i)  + 

6(ti)o  +  / 

;jj(*)i  =  +  B((x)t(u)o  +  (x)o('‘)i)  +  6<«)iK39) 

Note  that  the  rule  (22)  is  used  to  approximate  quantities  like 
(xu)0  and  (xu)(.  For  open- loop  operation,  (n)0  and  (u)i  are 
constant,  and  so  the  model  (39)  is  linear  and  time-invariant 
in  this  case.  A  simulation  of  an  open-loop  transient  is  shown 
in  Figure  8.  The  figure  compares  the  inductor  current  wave¬ 
form  of  the  actual  circuit,  the  inductor  current  waveform  of 
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Figure  8:  Comparison  o(  Refined  State-Space  Averaged 
Model  with  Traditional  State-Space  Averaged  Model 

the  state-space  averaged  model,  (i)o  which  is  the  dc  compo¬ 
nent  of  the  inductor  current  of  the  refined  model  (39),  and 
an  inductor  current  waveform  reconstructed  from  the  refined 
averaged  model  determined  by 

«(<)  =  (.•>-, +  (i)o(0  +  {i)i(*)e***.  (40) 


in  analyzing  resonant-type  converters.  Further,  the  approach 
offers  refinements  to  the  theory  of  state-space  averaging  that 
permit  ripple  estimates  and  offers  a  framework  for  analysis 
and  design  when  “small  ripple’  conditions  do  not  hold.  The 
method  may  find  applications  in  simulation  as  well  as  design 
since  it  is  considerably  easier  to  simulate  an  averaged  model 
than  a  switched  model. 


A  Time- Varying  Frequency 

As  briefly  discussed  in  Section  3,  the  analysis  based  on  (17) 
and  (18)  is  valid  if  w,  is  constant.  In  the  case  where  the 
frequency  w,(t)  is  time  varying,  it  is  appropriate  to  consider 
an  instantaneous  phase  function  defined  via 

!*(<)  =«,(<).  (41) 

In  this  framework,  the  Fourier  analysis  should  be  performed 
with  respect  to  the  basis  functions  {e;*4,<)}  for  0(t)  in  some 
interval  (0|  —  2*,0\]  rather  than  the  basis  The  time 

interval  for  the  analysis  now  depends  on  the  instantaneous 
phase  function  in  the  following  way.  Define  7(0  to  be  the 
duration  of  time  so  that  the  phase  function  varies  by  exactly 
2t  on  the  interval  (1  —  7(0<t]-  That  is,  for  a  particular  time 
t\,  7(0)  is  defined  by 


The  switching  frequency  has  been  selected  to  be  relatively 
low  so  that  the  ripple  is  emphasized.  As  can  be  seen  from  the 
figure,  tbe  additional  term  in  (39)  significantly  decreases  the 
error  l>e:»een  the  waveform  of  the  averaged  model  and  the 
underlyir.g  system.  Evidently,  it  is  possible  to  include  more 
terms  in  the  Fourier  series  representation  to  further  decrease 
approximation  error. 

The  refined  averaged  model  is  especially  useful  for  analysis 
and  control  design  in  the  case  where  ripple  is  not  small.  For 
instance,  the  eigenvalues  of  an  equivalent  continuous  time 
model  derived  from  an  exact  sampled-data  model  for  the 
present  example  can  be  determined  to  be 

Ai.j  =  jlog(\t.2{*(T)))  =  (—1.029  ±>1.328)  *10*. 

Note  that  4>(7)  is  the  transition  matrix  for  (38).  In  compar¬ 
ison,  the  eigenvalues  of  the  state-space  averaged  model  are 
given  by 

A",-*"  =  (-1.029  ±>1.230)  •  104 
while  those  of  the  refined  model  are  given  by 

Ar,/.w  _  (_i.029  ±>1.325)  *  104, 

(-1.029  ±>5.149)  *  104, 

(-1.029  ±>7.707)  *  104. 

One  pair  of  the  eigenvalues  of  the  refined  model  is  seen  to 
approximate  the  eigenvalues  of  the  underlying  system  much 
more  closely  than  those  of  the  usual  state-space  averaged 
model.  • 


=  *(M  (42) 

ff,-2»  =  *{t,-7(M}-  (43) 

Now  for  the  interval  (1  —  7'(t),  ij,  the  natural  generalization 

of  (17)  is 

*(f  -  7(0  +  s)  =  $:(x)*(«)«riM<'-"r'+*»  (44) 

for  s  6  (0,7(0).  In  this  case,  the  Fourier  coefficients  (x)*(f) 
need  to  be  defined  via 

(*>*(0  =  H  *{*(*))*■■’“  <# 

=  f‘  x(s)e--»w<*>u»(s)ds 

Jt-m 

rTi*) 

=  /  x(t  -  r(0  +  * 

Jo 

c j(t  —  T(t)  -f  s)da  (45) 

where  s(8)  in  the  first  line  of  (45)  represents  a  local  inverse 
map  from  phase  to  time.  An  easy  calculation  performed  after 
changing  the  variable  of  integration  in  (45)  by  a  =  7(f)  -  s 
yields 

£(*>*(0  =  ,{t  -  7(0)e-^*-«‘>U{f  -  7(0)7(0+ 

-  <7)e-'*'<'-#>a:,(<  -  <r)+ 

{w,((  -  <r)  -  >AvJ(f  -  -  a)c~'u,,~r)}do. 

(46) 

This  formula  is  equivalent  to 


0  =  -  7(0)e-^(,-r,,))w.{<  -  7(0)7(0  + 

(>M0  +  ((-->^.)x)t(0.  (47) 


5  Conclusions 

A  new  approach  to  averaging  in  power  electronic  circuits  has 
been  introduced,  and  the  method  has  been  seen  to  be  useful 
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Another  simple  calculation  yields 


^(0=1- 


«.(0 


(48) 


It  ><  now  clear  that  for  the  case  where  ui,  is  constant,  (47)  re¬ 
duces  to  (20).  Furthermore,  for  slowly  varying  w,(<),  the  last 
term  in  (47)  is  well  approximated  by  -jk ui,(r)t.  The  term 
x{t  -  r(/)}e-,M('‘r(0,w,{t  —  T(t)}T(t)  is  also  guaranteed  to 
be  small  if  w,(<)  is  slowly  varying  since  T(t)  is  given  by  (48). 

Another  important  consequence  of  the  formula  (47)  is  that 
there  is  no  direct  feedthrough  term  into  (*)*(!)  from  w,(l)- 
Hence,  in  the  case  where  the  frequency  ut,(t)  undergoes  a 
step  change,  we  observe  no  step  (or  impulsive)  behavior  in 
(x)t(i).  This  is  of  interest  in  the  examples  based  on  resonant 
converters. 
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